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Abstract 

The temperature dependence of the Casimir force between a real metallic plate and a metallic 
sphere is analyzed on the basis of optical data concerning the dispersion relation of metals such as 
gold and copper. Realistic permittivities imply, together with basic thermodynamic considerations, 
that the transverse electric zero mode does not contribute. This results in observable differences 
with the conventional prediction, which does not take this physical requirement into account. The 
results are shown to be consistent with the third law of thermodynamics, as well as being not 
inconsistent with current experiments. However, the predicted temperature dependence should 
be detectable in future experiments. The inadequacies of approaches based on ad hoc assump- 
tions, such as the plasma dispersion relation and the use of surface impedance without transverse 
momentum dependence, are discussed. 
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I. INTRODUCTION 



There are many corrections that one in principle has to take into account when calculating 
the Casimir force between two bodies; the corrections may come from finite temperatures, 
finite extensions of the plates used in the experiments, corrugations on the plates, etc. 
(Recent reviews of the Casimir effect can be found in Refs. 0,05 HI-) The correction that 
we will be concerned with in the present paper is the one coming from finite temperatures. 
For the most part, we will consider the temperature dependent Casimir force between a 
compact sphere of radius R and a plane substrate. The sphere is situated at a fixed distance 
a from the plane (a denotes the minimum distance between the surfaces). The sphere and 
the substrate are assumed nonmagnetic, but we consider the case where they may be made 
from different materials. We will moreover assume the proximity force theorem (41 to hold; 
this means that a must be much less than R. (For corrections to this see Refs. 0,1a Q-) O n 
experimental grounds it is evidently desirable to calculate the Casimir forces in a realistic 
way. We will here take advantage of the excellent numerical dispersive data to which we have 
access for the materials gold and copper (and also aluminum) (courtesy of Astrid Lambrecht 
and Serge Reynaud). We know how the permittivity e(*C) varies with imaginary frequency 
( over seven decades, ( G [10 11 , 10 18 ] rad/s. We use these data to calculate the forces at 
two different temperatures, namely at room temperature, T = 300 K, and at T = 1 K. 
The latter temperature is conveniently attainable numerically, and it can for all practical 
purposes be identified with zero temperature. (The lowest temperature that we actually 
tested was T = 0.2 K. If T becomes lower, we leave the frequency domain for our numerical 
dispersion data. It turns out that there are very small deviations between calculated values 
of the force for T = 1 K and for T = 0.2 K.) We obtain in this way a realistic picture of the 
finite temperature correction for these materials. 

It ought to be emphasized that we are not adopting the so-called modified ideal metal 
(MIM) model, which assumes unit reflection coefficients for all but the transverse electric 
(TE) zero mode [see Eq. ()3.2j) below]; rather, we are using real data together with the 
assertion (based on thermodynamical and electrodynamical arguments) that the TE zero 
mode is absent, as that is an isolated point which cannot be extracted from data alone. Our 
approach is that which we have followed in other recent papers [H, H, ©]. The absence of 
a TE zero mode contribution to the Casimir effect for a real metal was discussed in detail 
in Ref. and also in Ref. The result of this assumption, for instance, in the MIM 



model is the presence of a linear temperature term in the expression for the Casimir force 
between two planes, in the limit where aT <C 1. The calculation for a real metal yields a 
linear temperature correction for low, but not too low temperatures, so that for very low 
temperatures the force and the free energy have zero slope. By contrast, in the conventional 
(old) model for an ideal metal (IM) the TE zero mode is included, and it implies that this 
linear temperature term is omitted. We ought to stress here that at T = the mentioned 
difference between a MIM and an IM model goes away, as the contributions from the zero 
frequency TE mode as well as from the zero frequency TM mode become buried in an 
integral over imaginary frequencies from zero to infinity. 

The experiments of immediate interest for the present theory are the atomic force mi- 
croscopy (AFM) tests, performed in particular by Mohideen et al. [11]. A point that we 
ought to emphasize here is that previous analyses have most likely overestimated the accu- 
racy of the AFM experiments. Thus the recent paper of Chen et al. [l^, which is based 
upon a reanalysis of the experiment of Harris et al. (listed in 11]), claims an over- all exper- 
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imental precision to be at the 1 percent level. In that apparatus a gold-coated polystyrene 
sphere mounted on a cantilever of an AFM was brought close to a metallic surface and the 
deflection was measured. However, as discussed in Refs. jH [l3| , at the very short distance 
of 62 nm (the minimum distance) the force at a = 62 nm + 5 differs from the force at a = 62 
nm by more than 3.5 pN (the experimental uncertainty claimed by the authors) when 5 is 
larger than a few angstroms. This means that a should have been measured with atomic 
precision in order to correspond to the accuracy claimed. As for temperature corrections, 
these were found in [12] to be negligible, this being related to their acceptance of the plasma 
dispersion relation for the material. 

As the absence of the zero frequency TE mode has been controversial, we give in Sec. IIVI 
a discussion of this in view of recent work by Bezerra et al. [lij They argue that this mode 
should be present. In Sec.|V|we give additional support to our arguments by showing that for 
a pair of anisotropic polarizable particles the Casimir force can vanish in certain directions 
as the temperature increases towards oo, and although there are regions of negative entropy 
connected with the Casimir effect, there is no indication that thermodynamics is violated. 
Violation of thermodynamics is used as an argument by Bezerra et al. [lij to require the 
presence of the TE zero mode. The inadequacy of using a surface impedance approach 
without including transverse momentum dependence is briefly reviewed in Sec. IVII 

In this paper we put h = c = ks = 1. 



II. GENERAL FORMALISM AND DISPERSIVE PROPERTIES 

Let the sphere of radius R be nonmagnetic, and have a permittivity E\. As mentioned, the 
sphere is situated a distance a above a plane substrate; we let the nonmagnetic substrate have 
permittivity e-i. According to the proximity force theorem 0| the attractive force T between 
sphere and plane at temperature T can in the limit a/R <C 1 be given approximatively as 
the circumference of the sphere times the surface free energy density F in the parallel-plate 
configuration: T = 2irRF(a). Following essentially the notation of Ref. [|[ we can then 
write the force as 

T =lkY. \ ydy Wl-A™A™ C -*)+ln(l-A™A™e-^)]. (2.1) 

"° m =0 >^ m 7 

Here y is the dimensionless quantity y = qa, and 

Q=\/ki + (L ( m = 2nm/(3, 1 = 2na/(3, (2.2) 

k± being the component of k parallel to the plates in the parallel-plate configuration. Fur- 
ther, ( m with f3 = 1/T are the Matsubara frequencies, and 7 is the dimensionless temper- 
ature. Superscripts TM and TE in Eq. ()2.1j) refer to the transverse magnetic and electric 
modes. The prime on the sum means that the zero mode has to be counted with half weight. 
With the conventional Lifshitz variables defined as 

s = Ve-l+p 2 , p = q/( m , (2.3) 
we define the two kinds of A's, the reflection coefficients for a single interface, as 

A™ = £P ~ S , A TE = ^ (2.4) 
ep + s s + p 
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for each medium 1 and 2, respectively. If the two media are equal, Ai = A2 for each kind 
of mode, then 

(A™) 2 = A m , (A TE ) 2 = 5 m , (2.5) 

where A m , B m are the TM, TE coefficients defined in Ref. 0. Note that Si = 
y/ei — 1 + p 2 , S2 = — 1 + P 2 , with p = q/Cm being the same quantity in the two cases. 
The permittivities e(i( m ) are functions of the imaginary Matsubara frequencies ( m . In the 
general case where the media are dispersive, A™ and A TE depend both on p and on the 
Matsubara integer m. If the media are nondispersive, A™ and A TE are functions of p only, 
independent of m. 

As a general warning, we mention that the proximity theorem assumed here requires a/R 
to be very small. Thus, within the framework of the optical path method recently considered 
by Jaffe and Scardicchio j5j, believed to be more robust than the proximity approximation, 
disagreement with the latter approximation was found already when a/R became larger 
than a few percent, whereas the method they propose agrees accurately with the recent 
exact numerical result of Gies et al. ]m. 



A. Dispersive properties 

As mentioned in the Introduction, we will use accurate numerical data for the variation 
of e with frequency for two different substances: gold and copper. These data refer to room 
temperature measurements. For gold, the data are shown graphically in Refs. 0, El- For 
frequencies up to about 1.5 x 10 15 rad/s (note that 1 eV=1.519 x 10 15 rad/s), the data are 
nicely reproduced by the Drude dispersion relation 

^ = 1 + CKTW' < 2 ' 6 > 

where for gold the plasma frequency is uj p = 9.0 eV and the relaxation frequency v = 35 
meV. For ( > 2 x 10 15 rad/s the Drude curve however lies below the experimental curve. 

All the dispersive data, of which we are aware, refer to room temperature. Now, as 
we will be interested in the Casimir force also at low temperatures, we are faced with the 
problem of how to estimate the permittivity e(i(,T) under such circumstances. Numerical 
trials indicate rather generally that the Casimir force is very robust under variations in the 
input values for the permittivity, but at least this issue is a matter of principle. One possible 
way to proceed is to write the permittivity as 

and make use of the Bloch-Griineisen formula for the temperature dependence of the elec- 
trical resistivity p. This problem was discussed in Appendix D of Ref. 0. One may in this 
way estimate the temperature relaxation frequency to be, in eV, 

„(T) = 0.0847 g)/ JLiiL, (2.8) 

where O = 175 K for gold. This formula implies that e(i(,T) is somewhat higher for low T 
than for T = 300 K (cf. Fig. 1 in [8]), when the frequencies are less than about 10 14 rad/s. 
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It is instructive to compare with the parallel-plate configuration, where it is known that 
the most important frequencies for the Casimir force are lying in the region £a ~ 1. For 
a = 1 /jin it corresponds to ( ~ 2.5 x 10 14 rad/s. For smaller gap widths — where the metallic 
properties of the medium fade away and its plasma properties become more dominant — it 
follows that the most important frequencies become higher. Taking all things together, we 
expect that the influence from the temperature variation in v{T) is rather small. Some 
numerical trials that we have done support this expectation. 

Another important point to be mentioned here is that the Bloch-Griineisen argument 
sketched above neglects the effect from impurities. These give rise to a nonzero resistivity 
at zero temperature ^6[. This fact strengthens our argument for setting the contribution 
from the Casimir force from the TE zero mode for a metal equal to zero. The issue has been 
considered in detail also by Sernelius and Bostrom 10, 17]. As one can see from Fig. 2 in [10], 
there exists a temperature dependent contribution to the Casimir force from the TE modes. 
This temperature dependence occurs for low frequencies at low temperature and extends to 
higher frequencies with increasing temperature. However, the temperature influence fades 
away before the the first nonzero Matsubara frequency is reached. It is therefore permissible 
to neglect the temperature dependence in z/(T). What remains important in v is the constant 
term v{T = 0) ^ that is due to elastic scattering. The consequence is that 

C 2 K<)-1]^0, as C^O, (2.9) 

so from ()2.4j) A TE vanishes at m = 0. If one neglects this crucial constant v{T = 0), one can 
end up with a violation of the Nernst heat theorem jT^|. See also Bostrom and Sernelius 
for related discussion of these points. 
Therefore, in the following we will restrict ourselves to using the room-temperature values 
for v throughout, even when calculating J-(T) at different temperatures. Even with this 
simplification, we point out that the temperature dependence in the Casimir force turns up 
in a rather complex way. Namely, the temperature occurs at three different places: 

(i) in the prefactor in Eq. (|2.1j) ; 

(ii) in the lower limit of the integral; 

(iii) in the dependence of A™> TE and A™' TE on T via the Matsubara frequencies in the 
permittivity: e = e(i27miT). 

III. NUMERICAL RESULTS 

Figure Q shows, for a gold sphere and a gold plate, how the attractive force J 7 (a) varies 
with a in the interval from about 150 nm to 1 /xm, when R = 296 /im and T = 300 K. 
As mentioned, the empirical data for e(i() are directly usable as input in Eq. 1)2.1)1 . When 
a = 200 nm, the force is calculated to be 67.22 pN. 

A similar calculation can be made for a very low temperature, in order to show the 
magnitude of the temperature influence in the force. We assume throughout, as mentioned 
earlier, that v = 35 meV. Using MATLAB we found T = 1 K to be a numerically stable 
and reasonable lower limit. This temperature is moreover low enough to be identifiable with 
T = for all practical purposes. 

We ought to stress that we choose to perform the T m calculation numerically, inserting 
realistic data for e(*C). This is in principle different from the conventional T = calculation 
for an idealized metal, where one simply puts e = oo for all frequencies. (Recall that at 
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FIG. 1: Force T between a gold sphere and a gold plate versus gap a, when T = 300 K and R = 296 
jum. 

T = the difference between a MIM and an IM model goes away, because of the very close 
spacing between the Matsubara frequencies.) 

Rather than showing the calculated result for J 7 (IK) explicitly, we show in Fig. [21 the 
difference between the forces, AJF, defined as 

A:F = :F(300K) -:F(1K) = 1^(1 K)| - 1^(300 K)|. (3.1) 

An important property seen from this curve is that AJ 7 is positive. The force is thus weaker 
at room temperature than at T = 0. This is the same effect as was found in Fig. 5 in 0. 
This behavior is thus a consequence of Lifshitz' formula plus realistic input data for the 
permittivity; there are no further assumptions involved. When a = 200 nm, we find AjF= 
2.54 pN, which means that the force is reduced by 3.6 percent compared to the T = case. 

For larger distances, a = 400 nm, the temperature effect becomes larger. Thus for a = 400 
nm the force is 9.38 pN at T = 300 K and 10.19 pN at T = 1 K, yielding a 7.9 percent 
reduction at room temperature. At a = 1 /im, the corresponding numbers are 0.59 pN at 
T = 300 K and 0.73 pN at T = 1 K, which means a 19 percent reduction. This agrees in 
magnitude well with the temperature corrections in the case of parallel-plate geometry, as 
is seen from Fig. 5 in Ref. 0. 

Admitting an error of 10~ 8 in the m summation in Eq. (J2.1)) . we found the necessary 
number of terms to be in excess of 34000 in the case of the lowest separation investigated 
numerically, a = 50 nm (not shown in the figure). When a = 200 nm, about 11000 terms 
were required. At larger separations the necessary number of terms became considerably 
reduced; thus the case a = 1 fim corresponded to about 2700 terms. 
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FIG. 2: Force difference AJ 7 = |.F(1K)| — |.F(300K)| between a gold sphere and a gold plate, 
versus gap a for R = 296 fim.. 

The calculation of the force between a gold sphere and a copper plate gave very similar 
results. Thus for a = 200 nm the force was 67.19 pN at T = 300 K and 69.75 pN at T = 1 K, 
corresponding to a reduction of 3.7 percent at room temperature. At a = 1 /im, the forces 
turned out to be the same (to the accuracy of two decimals) as in the Au-Au case. 

In Fig. El we show, by the full line, how the force between a gold sphere and a copper plate 
varies with a, at T = 300 K. The curve is calculated from Eq. (|2.1j) . using the empirical data 
for these two materials directly. Our reason for giving this curve anew, in spite of its identity 
with the curve in Fig. ^ for all practical purposes, is that we have supplied the following 
new element, namely at the bottom of the figure we show explicitly the contribution from 
the m = term to the force. Namely, for the MIM model in which Aq — 1, Bq — 0, A m = 
B m = 1 for m > 1 [H, we have for the free energy 
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showing the m = contribution in the first term 
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meaning that the corresponding force contribution 2irRF 
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FIG. 3: Force T between a gold sphere and a copper plate versus a, when T = 300 K and R = 296 
/im. Bottom curve is the m = contribution. 



From the figure this term contributes an increasing part of the total force for increasing a, 
and gives the full contribution in the classical limit a — > oo. 

We would like to make a direct comparison of our prediction with the recent experiment 



of Decca et al. |20l. l21l|. This experiment measured the Casimir force between a gold sphere 
and a copper plate by means of a microelectromechanical torsional oscillator, for separations 
in the range 0.2-2 /zm. The radius of the sphere was 296 ± 2 //m, thus the same radius 
as we have assumed above. Their measured values are shown, for instance, in Fig. 3 in 
Ref . • However, the scale of their figure for the total force prevents comparison with our 
numerical results with sufficient accuracy to draw firm conclusions about the magnitude of 
the temperature dependence. 

But Fig. El in Ref. [icj also gives a comparison with theoretical values by plotting the 
difference. These theoretical values are evaluated at T = 0, but as we do, they also use a 
Drude model to obtain the dielectric constant in the limit of low frequencies. Then they 
find a small temperature correction with the same sign as we have theoretically predicted. 

From Fig. |3] in Ref. we estimate AT to be around 1 pN for a = 200 nm, and to 
be too small to be discernible for large gaps, a > 600 nm (the scatter of the experimental 
points is considerable). This can be compared with our Fig. |21 where we calculated AT to 
be 2.56 pN for a = 200 nm and 0.14 pN for a = 1 /im. This reasonably good agreement 
between experimental and theoretical results is encouraging, and it indicates that our finite 
temperature calculations are on the right track. 

We should also mention that Ref. |2l| refers to dynamical measurements that are claimed 
to rule out our results. However, we believe that the systematic theoretical and experimental 
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uncertainties in this measurement are larger than estimated by those authors. In particular, 
we re-emphasize the uncertainty of measuring the Casimir force due to uncertainties in 
estimating a systematic shift of position as earlier discussed by Iannuzzi et al. [13| arid 
Milton as mentioned in the Introduction. 



IV. COMMENT UPON THE TRANSVERSE ELECTRIC ZERO MODE 

Bezerra et al. 3 claim that the Drude dielectric function for metals cannot be used in 



the theory for the Casimir force. Their opinion is that it violates the Nernst theorem in 
thermodynamics and is furthermore ruled out by a recent experiment. Instead the plasma 
relation should be used for the dielectric function. The latter implies the presence of a 
transverse electric mode at zero frequency besides the static dipole-dipole interaction; the 
latter being the zero-frequency limit of the transverse magnetic mode. Including such an 
transverse electric mode adds a term linear in temperature to the Casimir force by which it 
is increased by a factor of two in the high temperature limit. 

These authors point to the standard theory, sketched above in Sec. Ill Al where the relax- 
ation parameter v varies with temperature T and vanishes at T = 0. This is then a situation 
where in principle a transverse electric zero mode might be present in the Casimir force 
although it should not be present according to Maxwell's equations of electromagnetism. 
However, its presence for v = 0, but not otherwise, would make physical phenomena dis- 
continuous as then v = will give something different from taking the limit v — > 0. Also 
we question whether a possible vanishing of the relaxation parameter at T = can dictate 
the behavior for T > where v > 0. (As a remark we here note that strictly speaking 
the statistical mechanical derivation is exact only for dielectric functions independent of 
temperature, i.e., when induced dipole moments are harmonic oscillators. A T-dependence 
reflects anharmonicity. But we do not expect that this is of crucial significance here.) 

As a result, Bezerra et al. [14] conclude that the Drude dielectric function is thermody- 
namically inconsistent and cannot be used to calculate the thermal Casimir force for metals. 
However, we on the contrary have shown that it is thermodynamically consistent leading to 
zero entropy at T = in accordance with the Nernst theorem (if. Thus, for T > a negative 
entropy contribution due to electromagnetic interaction between media is allowed as long as 
the total entropy is positive. We have earlier studied a simplified model in this regard 
The model consists of three harmonic oscillators. Two of them are the analogues of the two 
media that interact via the electromagnetic field represented by the third oscillator. 

In the beginning of their Sec. Ill, Bezerra et al. 3 state that our derivations were based 



upon a constant relaxation parameter v. However, we did not require such a limitation, 
only that z/(0) is finite. And as the authors further note, the value of v is commonly very 
small, but nevertheless finite, at T = due to impurities, and as a result the entropy 
becomes zero at T = as required. Then they write that the negative entropy we found 
at T = 5 x 10~ 4 K is a violation of the Nernst theorem. But such a negative perturbing 
entropy is not a violation of the Nernst theorem since T is finite as already discussed above. 
The further discussion about relaxation time, impurities, and surface impedance does not 
invalidate the use of a finite ^(0). 

These authors further remark that we consider nonzero wave vectors for which the 
reflection coefficient r\(0,kj_) = 0. But this does not mean that reflection properties are 
different for the fluctuation field compared to real photons as we deal with only one such 
quantity; and as discussed above, thermodynamics is not violated. (See Sec. |V3 below for 
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further discussion of dependence.) 

They further note that a material with dielectric constant e = 100, like a metal, gives a 
Casimir force that decreases with temperature in some interval and thus implies a negative 
entropy contribution. They conclude that such a material cannot exist as it would violate 
thermodynamics. Then they claim that real media with such large e are commonly polar for 
which e rapidly reduces to its optical value connected to its electronic polarizability. So with 
plate (or plate-sphere) separation of 1 fim the Casimir force again becomes monotonically 
increasing with T. However, this does not preserve the monotonic character of the Casimir 
force in general because one can just increase the plate separation. Non- monotonic behavior 
will then reappear when this separation exceeds a wavelength corresponding to the relaxation 
frequency (where e decreases rapidly). 

In Sec. IV of their paper, they again conclude that the use of the Drude dielectric function 
in the Lifshitz formula violates the third law of thermodynamics (the Nernst heat theorem). 
This conclusion, stated as a rigorous proof, is made on the basis that the relaxation param- 
eter will be much less than the Matsubara frequencies. But this is not a rigorous proof as 
the relaxation parameter will violate this inequality sufficiently close to T = (assuming 
v(0) > 0). Furthermore the Drude dielectric function does not predict a linear temperature 
correction to the Casimir force all the way down to zero temperature. As we have shown 
earlier the Casimir force flattens out and becomes independent of temperature at T = in 
accordance with thermodynamics [H0]. (But the sharpness of this flattening increases with 
decreasing v.) 

We further remark that the use of the Drude dielectric function is consistent with the 
e — > oo limit (for large separation of the plates). Use of the plasma model {y — 0), however, 
will give a discontinuous jump of the force in this limit. Such a jump is not expected for 
a continuous change in a physical parameter. Also for real metals the e effectively will be 
finite due to the finite size of the plane-sphere configuration. (Here one can note that a 
medium consisting of separate metal spheres of finite size will be like a polarizable medium 
with finite polarizability and dielectric constant.) 

Finally the authors in Ref. 0] conclude that use of the Drude dielectric function ()2.6j) 
is in contradiction with experiment. However, so far as we can see, experiments performed 
at a single temperature are at present too uncertain to draw conclusions about temperature 
variations. This seems even more evident from Ref. 2 • There detailed analysis of exper- 



imental uncertainties are performed and various corrections for very short separations are 
estimated. They find an uncertainty of 1.75 % at 95 % confidence level for 62 nm separation. 
This uncertainty increases to 37.3% for 200 nm. Earlier it has been remarked by others 
[l3l | that such experiments are very sensitive to accurate determinations of plate separation 
since the plate-sphere Casimir force is proportional to the inverse cube of separation dis- 
tance. However, the authors assert that they avoid this problem by making a least squares fit 
of the resulting data by which zero separation is pinpointed with an uncertainty of 0.15 nm. 
But in view of the uncertainties of the experiment this does not seem to resolve the disputed 
temperature dependence as heavy weight is put on the shortest separation of 62 nm where 
the force is by far the largest and changes most rapidly. Thus high apparent precision can 
be obtained for this separation (1.75%) while for larger separations the uncertainties are 
rapidly increasing until they become larger than the magnitude of the disputed thermal 
effect. According to the authors of Ref. l3| the thermal effect in dispute is about 1-2% for 
62 nm. 

Very recently, a paper has appeared 2^| giving the microscopic theory of the Casimir 
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effect. These authors convincingly demonstrate that the TE zero mode cannot contribute, 
although the plasma model gives such a contribution. 

V. ANISOTROPIC PARTICLES WITH NEGATIVE CASIMIR ENTROPY 

As mentioned above the appearance of a negative Casimir entropy for metals in a region of 
nonzero temperatures has been disputed with the claim that it violates the Nernst theorem 
of thermodynamics 0|. However, as we argue this negative entropy region does not imply 
violation of thermodynamics since the Casimir free energy is a perturbing one and is not the 
total free energy of two interacting systems. Many such examples are known in statistical 
mechanics, including that of three interacting oscillators discussed in Ref. ||. To illustrate 
this point, we here will consider a pair of particles with strong anisotropic polarizability and 
thus polarizable only in the ^-direction, e.g., they may be metal needles. The result for a 
pair of particles with isotropic polarizability is well known and was rederived in a novel way 
by Brevik and H0ye j23|. 

The latter derivation is easily generalized to anisotropic particles. Using the path integral 
formalism the dipole-dipole interaction of the radiation field is given by Eq. (15.9). (Here 
and below the numeral I refers the the equations of the reference mentioned.) Equation 

(15.9) gives the interaction energy 

oo 

0S(12) = V [^ n (r)D Cn (12)+i; ACn (r)A Cn (12)]a 1Cn a* Cn (5.1) 

n=— oo 

with 

D Cn (12) = 3(r • a 1? J(r ■ a^J - a 1?n • a^ n , (5.2) 

A Cn (12) = a lc „-a; Cn (5.3) 

where the Matsubara frequency ( n = 2im/f3. The dipole radiating fields ip are given in 

(15.10) . Hats denote unit vectors. The a^ n are unit vectors of the Fourier-transformed fluc- 
tuating dipole moments along the "polymer" path and r] — > is the discretized step length 
along the "polymers" that represent quantized particles in the path integral formalism. 

The Casimir free energy times —2(3 is now obtained from the average of expression (|5.1|) 
squared in accordance with (13.5). For the isotropic case one has the average ^(a^a^) = 
2>aQ n where a^ n is the polarizability. This follows from Eq. (15.3). For the strongly anisotropic 
situation to be considered below one likewise will find (z denotes the ^-component) 

as the only nonzero average as polarizations in x- and ^/-directions are zero in the present 
case. 

Furthermore let the positions of the particles relative to each other be such that the 
z-component of f equals 1/VS. With this the D^ n (12) will vanish as polarizations are 
present only in the z-direction. (With this relative position the corresponding electric field 
is transverse to the z-direction and thus to each of the dipole moments, i.e., there is no 
interaction connected to the D^ n (12)-term.) Thus only the A^ n (12)-term remains. So with 
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polarizations restricted to the z-direction Eq. (15.14) turns into 

~P F =\ll<ASr) (5.5) 

n 

or with (15.10) inserted the free energy is 



where 




as given by (15.12). 

As usual Eq. (|5.fij) gives a negative free energy. However, in the classical limit T — > oo the 
Casimir free energy (|5.6|) . and thus the corresponding force, both vanish since only the n = 
term will contribute. Furthermore with this free energy the contribution to the entropy must 
be negative (or at least mainly negative) as S = —OF/ dT, because F must have a generally 
positive slope, but S = at T = as it should. 



VI. SURFACE IMPEDANCE 



Most recently, Mostepanenko et al. |14L |24 , apparently conceding that their arguments 



favoring the plasma model over the Drude model for the dispersion relation for real metals 
could not be supported either thermodynamically, electrodynamically, or experimentally, 
have asserted that for real metals one should use in the reflection coefficients in the Lifshitz 
formula not the bulk dielectric permittivity but the surface impedance. Indeed there is much 
to be said for using the latter. However, in may be shown in general that there is in fact no 
difference between the reflection coefficients computed using either description There 
is a one-to-one correspondence between the permittivity e and the surface impedance Z, 
which is given by the ratio of the transverse electric and magnetic fields at the surface. This 
correspondence, however, necessitates in general that both quantities possess dependence 
on the transverse momentum kj_. As optical data strongly suggest that this dependence is 
usually negligible in the permittivity, a strong dependence on k^ is required in Z H 

zTE(c '^ ) = -^r4Y=w- (e- 1 ) 

The vanishing of 1 + Zq/( at ( = demonstrates again that the TE zero mode does not 
contribute to the Casimir force. (Note that this vanishing happens in the Drude, but not the 
plasma model.) In contrast, Ref. 0, |24| completely disregard this transverse momentum 
dependence and moreover make an ad hoc extrapolation from the infrared region to zero 
frequency |]| . The inadequacy of neglecting the transverse momentum dependence has been 
stressed by Esquivel and Svetovoy 25]. 
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VII. CONCLUSIONS 



In this paper we have sharpened our arguments in favor of using real data for the dielectric 
functions to apply the Lifshitz formula to calculate the force between metal surfaces, in 
particular between a spherical lens and a flat plate. In principle, one can also use the 
surface impedance to calculate this force, and although optical data is lacking for the latter, 
such use should yield the same result. In contrast, the procedures advocated in Refs. 0, 
lil . ill . I20I . I21L l24j contain ad hoc elements and assumptions. 

We show both by direct computation, and through analogous models, that there is no 
conflict with thermodynamical principles, in particular with the Nernst heat theorem (the 
third law of thermodynamics). Especially important is the demonstration that the entropy 
necessarily vanishes at zero temperature. Claims that experimental limits on Casimir forces 



preclude our temperature dependence |20j, |21[ are, in our opinion, not justified, since the 



accuracy of the current experiments does not match their precision, especially due to the 
impossibility of determining the shortest separation distances accurately. Undoubtedly, it 
will take dedicated experiments involving different temperatures to reveal the true tempera- 
ture dependence of the Casimir effect. One idea might be to measure the difference between 
the Casimir forces for the same value of a at two different temperatures, for instance 300 
K and 350 K. Such a difference is directly measurable, in principle. To our knowledge this 
idea was originally proposed by Chen et al. |26|, and it was further discussed in Ref. 
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